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In this paper we take up the quantal two-centre problem where the Coulomb centres have arbitrary 
positive charges. In analogy with the symmetric case, treated in the second paper of this series of 
papers, we use the knowledge on the quasiclassical dynamics to express the contour integrals in 
the first- and third-order approximations of the phase-integral quantization conditions, given in the 
first paper of this series of papers, in terms of complete elliptic integrals. For various values of the 
distance between these charges the accuracy of the formulas obtained is illustrated by comparison 
with available numerically exact results. 
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1. INTRODUCTION 



In this third paper in a series of papers concerning the two-centre Coulomb problem we take up the general case 
in which the positive charge numbers Z\ and Z^ of the two Coulomb centres are arbitrary. The study of such 
general systems is of considerable importance in the field of molecular physics and elementary particle physics. For 
example, the calculation of eigenvalues and electronic wave functions for one-electron diatomic molecules with fixed 
internuclear separation is the starting point for an accurate description of molecular vibrations and rotations and of 
ion-atom scattering jl], g Qj. Also the calculation of radiative transition probabilities for a ir~ meson moving in the 
Coulomb field of two fixed nuclei Q is a good example of the two-centre Coulomb problem dealt with in elementary 
particle physics. Such calculations are of physical interest in connection with experimental |6[ 0, ||] and theoretical||] 
research on the absorption of tt~ mesons stopped in substances containing hydrogen. 

In the second paper in this series |l(| the symmetric case, Z\ = Z2, was considered. Using for the two-centre 
Coulomb problem the general phase-integral quantization conditions derived in the first paper JTl| ], the relevant 
contour integrals for the first- and third-order approximations were expressed in terms of complete elliptic integrals 
so that numerical evaluation of energy eigenvalues and separation constants can easily be carried out. The evaluation 
of the various contour integrals was facilitated through suitable transformations of the £- and 77- variables, which can 
be related to the quasiclassical motion of a particle. In the symmetric case, Z\ = Z%, the square of the base function 
Q 2 (rj) has a symmetry, as discussed in jlo|]. As a result of this, the evaluation of the quantities a = (3, L, L' and 
K pertaining to the 77-quantization conditions were performed with the use of particular transformations in a rather 
simple way, while for the quantity L in the ^-quantization condition rather general transformations are necessary. 
However, in the general case, where Z\ may be different from Z2, no such symmetric structure exists for Q 2 (r]), and 
both Q 2 (ri) and Q 2 (£)j given by eqs.(3.2a,b) in jlj)] , can have, besides the poles, the following structure of the zeros: 

1. Case A = \m\ ^ 0: 

(a) Four real zeros 

(b) Two real and two complex conjugate zeros 

2. Case A = 0: 

(a) Two real zeros 

(b) Two complex conjugate zeros 
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As a consequence, the evaluation of the contour integrals in the cases of the ^-equation and the 77-equation are 
analogous, except that certain coefficients change in the different subcases. This allows one to use the "universal" 
functions ij( 2n+1 ) and H^ 2n+1 \ n = or 1, valid for the (2n + l)th-order approximation, which were defined in 
sec. 3.1.2 and sec. 3.2.2 of ||Io| . These functions depend on five parameters v\ , V2 , 1^3 , g and k, which are expressed 
in terms of the zeros of either Q 2 (7 7) or Q 2 (£), and on the parameter C and C in the base function or Q(rj), 

respectively; see eqs. (3.2a,b) in pM. Different choices of these parameters allow one to express the relevant quantities 
a, (3 (which may now be different from a), L,L',K and L in the first- and third-order phase integral approximation 
(apart possibly from a sign or a constant factor) as the appropriate "universal" function or its real or imaginary part 
with the appropriate parameters inserted. When the contour integrals are evaluated, one can solve the quantization 
conditions to obtain the energy levels accurately. 

In principle one can specialize the results of the general case when Z\ may be different from Z 2 to the particular 
case when Z\ = Z2 in order to obtain the results of Jl(|, but in practice this is cumbersome due to the different 
transformation formulas used for the 77-part of the quantization conditions in JlO[ . For the £-part such a specialization 
implies only that Z\ + Z2 is replaced by 2Z\ but no essential simplification. 

The plan of the present paper is as follows. In sec. 2 we express a,(3,L,L',K and L in terms of the "universal" 
functions appropriate for the subbarrier case with A = \m\ ^ 0. Then a similar treatment is given for the superbarrier 
case. In sec. 3 an analogous procedure is applied to the case A = 0. Finally, in sec. 4 a detailed numerical analysis 
of the phase-integral quantization conditions is carried out for Z\ = 1 and different values of Z2 (=2, 5 and 8), and 
comparision is made with existing numerically exact results for the energy and the reduced separation constant. 



2. CASE A = \rn\ / 

In this section we utilize for the case of four zeros of Q 2 (rj) or Q 2 (£) the "universal" functions H^ 2n+1 \ n = or 1, 
introduced in JTo| and related to the (2n + l)th-order contribution, and the "universal" functions H^ 2n+1 \ n = or 
1, which will be introduced in sec. 3.2.2. In the different cases one finds, as explained in the introduction, that the 
forms of the evaluated contour integrals are similar, except for changes of the parameters. One obtains the expressions 
for the "universal" functions H^> and by integrating one specific integral in each case explicitly and 

then modifying the definition of the parameters in these functions to obtain the other required quantities. 

2.1. Four real zeros of Q 2 (r]) or Q 2 (£) 
2.1.1. The quantities a, (5 and K pertaining to the n-equation: Subbarrier case [Fig. 3(a) in Ref. 11] 



Denoting the zeros 771 , 772 , ^3 and 774 by a, b, c and cZ, we write the base function for this case as 

[(77 - a)(b -r])(c- 77) (d - 77)] 

1—7/ 

Using the transformation on p. 103 in 113] , we obtain 



Q( V ) = p ^-^-^-w-'M' . (2.1) 



a — dvisn u . . 
h^i—' (2-2) 

the parameter v 2 being defined in (2.4) below. Noting that the loop a — > b —> a, that is 77! — > 772 — > 77! in the 77-plane, 
denoted by T a ^, which represents the contour A a in Fig. 3(a) of 
one finds that the first-order contribution to a is 

aW = i / Q( V )d V 



11|, corresponds to — > K — > 2K in the 77-plane, 



A,, 



r 



Q(r))dr) 
P f ZK ( df) \ 2 du 



2.9 Jo \du J 1 — 77 2 

2-p{a-d) 2 v\ f 2K sn 2 u(l - sn 2 u)(l - k 2 sn 2 u)du 



g(l — a 2 ) J (1 — v 2 sn 2 u) 2 (\ — v 2 sn 2 u)(\ — v 2 sn 2 u) ' 



(2.3) 
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where 



d-b' 



[(d-6)(c-o)]i 



(d-c){b-a) 
(d-b){c-a)' 



(2.4) 



(2.5) 



Decomposing the integrand in (2.3) into partial fractions and evaluating the integrals by means of recurrence 
formulas in Jl2j, we obtain 



*W =H^(v 1 ,u 2 ,u 3 ,g,k,C), 



(2.6) 



where is the "universal" function defin ed i n eqs.(3.16a,b), (3.17a-c) and (3.18) of jio) but with C replaced by C 
and the parameters vi, i = 1, 2, 3, given by (2.4), and g and k 2 given by (2.5). The evaluation of the quantity in 



the present case is thus similar to the evaluation of the quantity £w in sec. 3.1.2 of 1 1 
Proceeding in a similar way, the third-order contribution to a is found to be 



= - 
2 



-C 



l-rfj 2Q(ry)(l-?7 2 ) 
= H&\vi,v 2 >vz,g,k,C), 



1 

16 



'""if I "" 



(2.7) 



where is the "universal" function given in eqs. (3.16c,d), (3.19a-d) and (3.20a-d) of Q but with C replaced by 
— C and the parameters vi, i = 1, 2, 3, given by (|2.4|), and 5 and fc 2 given by (p^). 

In a similar manner as above we proceed for the calculation of the first- and third-order contributions to (3. The 
loop d — > c — ► d, that is 7/4 — > 7/3 — > 774 in the 77-plane, denoted by r^c and represented by the contour in Fig. 3(a) 
of |ll|| , corresponds to — > K — » 2if in the u-plane, and hence the first- and third-order contributions to (3 are 



(i) 



Q(r/)di] 



Q{vi)dri 
H^(i>i, v 2 ,v^g,k,C), 



and 



where now 



^ 3 ) = J ff( 3 )(r/ 1 , V2, v 3 ,9,k,C), 



and g and fc 2 are given by (2.5). 
Similarly one obtains 



where now 



d-b' 

ttKq 
■KK2 

c-b 
d-b' 



1 -c 



(^1,^2, v3,g,k,C), 
H ( - 3, >(vx,V2,vz,g,k,C), 



1 



1 



1-d 



and g is the same as given in (2.5) and k is now given by 



(c-b)(d-a) 



(2.8a) 
(2.8b) 

(2.9) 

(2.10a) 
(2.10b) 

(2.11) 
(2.12) 



(d-b){c-a) 

According to eq. (3.18a) in |ll|] the integrals a' and f3' for the contours A Q / and Ap> in Fig. 3(a) in |ll[ are obtained 

f and & = 13 + ?f 



from the formulas a 1 = a + 4f and j3' = (3 + A7r 
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2.1.2. The quantities L and V pertaining to the £- equation [Fig. 1 in Ref. 11] 

The formulas for the first- and third-order contributions to L and L' are derived and presented in subsection 3.f .2 
of ]To| , and they remain unchanged in the present case. 

2.2. Two real and two complex conjugate zeros of Q 2 {rf) or Q 2 (£) 

When there are two real and two complex conjugate zeros of Q 2 (rj), the situation of either Fig. 4(a) or Fig. 2 in 
pT| may occur. The latter situation has, however, so far not appeared in our applications, and therefore we disregard 
it in our treatment below of the 77-equation. 

2.2.1. The quantities a,f3, K, L and L' pertaining to the rj-equation: Superbarrier case [Fig. 4(a) ° r Fig. 2 in Ref. 11} 

Denoting the real zeros of Q 2 (rf) by 771 = a and 774 = d and the complex conjugate zeros 772 and 773 by c and c*, we 
have the base function 

1 — ?7 2 

Defining 

c = b%— iat, c*=b±+iax, (2-14) 

A=[(a-h) 2 + a 2 ]i, (2.15a) 

B = [{d-b l ) 2 +a\]^, (2.15b) 
and using the transformation on p. 133 in [Q, we get 

aB + dA + (dA - aB)cnu 
V = A + B + (A- B)cnu ' ( ' ' 

Here we exploit the fact that the Jacobian elliptic functions are doubly periodic, one of the periods being complex. 
Thus the loop d — > a\ + ib\ — > d, that is 774 — > rj 3 — > 774 in the 77-plane, denoted by T^c* , for the contour in Fig. 
4(a) of [[n], corresponds in the 7j-plane to the path — > K + iK' — > 2K + 2iK', where K and K' are complete elliptic 



integrals of the modulus k, given in ( 2T8e) below, and of the complementary modulus k' = yl — fc 2 , respectively. 
Making use of the transformation ( 2.1q ), we obtain for the integral in the first-order expression for [3 



\ f Q(v)dv = \f Q(r))dr) 



p(i>l — v 2 )(i>i — V3) r 2K+2lK ' sn 2 udn 2 udu 



2g J (1 + ^icn7i) 2 (l + V2cnu)(l + v^cnu) ' 



(2.17) 



where 



A — B , 



{l + d)A+(l + a)B> y ' 
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f 3 = 



(1 - d)A - (1 - a)B 
(l-^A + il-^B' 

1 
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<AB' 



(a-df-(A-Bf 



k = 



AAB 



(2.18c) 
(2.18d) 
(2.18e) 



Similarly one can treat the corresponding integrals in the first-order expressions for a, K and L. When one then 
evaluates the integrals containing the elliptic functions [cf. (2.17)], one finds that all these integrals can be expressed 
in terms of a "universal" function that is given by 



H^iui, f2,V3, g,k,C) = --{^2 

9 Uf 



2k 2 



n 



, k )+ {v\ - 2k 2 )K{k) 



-2v\E(k) + vxk(2j + 1)tt 



v\(\-2k 2 ) + 2k 2 v x 
[(1 - v\)(k 2 + iyfk' 2 )]i^ 



J7T 



i=l 



2k 1 



n 



,k'\ -2k 2 K{k') + 2E(k') 



where C%, C 2 and C3 are given by eqs. (3.37a-c) in flCfl , that is, 



Ci = 



Co 



(i>i - v 2 )(v x - v 3 ) 



{v\ - v z )v\ 



{V X - V 2 )(v 2 ~ ^3)' 



c. 



(1/1 - I/ 2 )ff 



(^1 - ^3)(^3 - V 2 ) 



(2.19) 

(2.20a) 
(2.20b) 
(2.20c) 



and 



1 



fc 2 (l - i^)Jf (fc) + v}E(k) - (k 2 + ufk 2 )II 
,,2 



itf - 1 



2kvf 
i 



k\»l - 1) - ^ ) (2j + 1)tt + fc[(l - v 2 ){k 2 + k 2 v 2 )\ijv 



+ 72 \-{k 2 + v 2 k 2 ) 



U (jh?' k ')- K{k,) 



-E(k') \,i = 1,2,3, 



(2.21) 



j being an integer —1, or +1, depending upon whether the quantity (3^, or oS l \ respectively, is evaluated, and 
originating from terms tan -1 (sdu) and cos -1 (dnu) while applying the limits of integration. Note that the quantity 
given by eq. (3.36a) in ]l0|] , is just —2ReH^ with j = 0. For if, which is expressed in terms of the imaginary 
part of fiW, the value of j, which appears only in the real part of H^ x \ does not matter. 
Similarly we also introduce the "universal" function given by 



H^{ VuV2 , V ^g,k,C) = -f- 

lop 



-4C* 



1 



^1(^2 - V3) 2 



v 2 v-i{v\ - v-s)(v 2 - V\) 



K(k) 



(v\ - vz)(v\ - v 2 ) \ 



(xK{k) + YE{k) + ^-K 



+1 



■ 9 



16p 



-AG 
1 



v\{v 2 - v z f 



K(k') 



(y\ - vz)(y\ - v 2 ) 



v 2 vz{v x - v z )(y 2 - v x ) 

{XK{k') +Y[K{k') - E(k')]} 



(2.22) 
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where C is the parameter in the square of the base function Q(rj) in eq. (3.2b) of pj, and 

l + 4k 2 3 + 4fc' 2 n n . fc' 2 . 4l2 . 2 

X = + {v 2 + 2 Vl v 2 +1v x v % + v 2 Vi)- gjs(17- 4fc 2 K^3 

-2^(^i + *a, + 1/3) + 2^ 3 + ^ 2 f — + ^ J , (2.23a) 

Y = -^L(l + 8fc 2 -8fc 4 ) + ^(2£ : 2 -l)( l / 1 2 + 2^ + 2^ + 1^3) 

+^P(17-8fc 2 + 8fc 4 ), (2.23b) 

Z = 2^1(1/1^2 + ^1^3 +4^2^3)5 (2.23c) 



where the parameters i^i , v 2 , v% , <? and k are defined in (2.1Sa-e) 



Then the first- and third-order contributions to the quantities a, (3 and K (= irK) are 

= Rei? (1) with j = 1, (2.24) 

= ReiJ (1) wiift j = -1, (2.25) 

tt^o = -2ImiJ (1) , (2.26) 

a (3) = Re# (3) with j = 1, (2.27) 

/? (3) = ReiJ (3) wiift j = -1, (2.28) 

7r^ 2 = -2±m7J (3) , (2.29) 



where v\, v 2 , V3, g and k 2 are still defined by (2.18a-e). Since j appears only in the real part of H^ 2n+1 \ we need not 
specify a value of j in (2.26) and (2.29). 

The integrals a' and (3' associated with the contours A a i and Ap> in Fig. 4(a) in (ll) are obtained from a and (3 by 
means of the relations (3.18a) in jyj, that is, a' = a + 4? and /?'=/?+ 4p. The integrals L and i' associated with 
the contours and Ar,> in Fig. 4(a) in |Tl"| ] can be obtained from a and 3 by means of the formulas L = a + f3 and 
L' = L + \m\. 

2.2.2. The quantities L and V pertaining to the ^-equation [Fig. 1(a) in Ref. 11] 

The formulas for the first- and third-order contributions to L and L' remain the same as the ones presented in 
subsection 3.2.2 of Q. 

3. CASE A = 

3.1. Two real zeros of Q 2 (v) or Q 2 (0 

One should be able to obtain the formulas pertaining to the case A = by considering the limits of the "universal" 
functions and when a — > — 1 and d — > +1 in the case of the ry-equation. These specialization procedures 
are, however, much more cumbersome than the direct calculation of the quantities in question, and so we do not carry 
them out here. Instead we shall evaluate these quantities directly, and therefore no "universal" functions will appear 
in subsection 3.1. 
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3.1.1. The quantities a, (5 and K pertaining to the rj-equation: Subbarrier case [Fig. 3(b) in Ref. 11] 

We denote the two real zeros of Q 2 (rj) by 772 = b and 773 = c and use transformations on p. 103, p. 120 and p. 
112 in for calculating the first- and third-order contributions to a, j3 and K, respectively. Here the base function 
reads as 



Q(v) =p 



(rj - b)(c- rj) 

1-T? 



The first- and third-order contributions to a are 



and 



where 



q(1) = ^^ E ( k ) ~ ^ + " 2 ) K ( k ) + ( fc2 - ^ 4 )n(^ 2 , fc)], 



(3) gCK(k) 
or ' = 



p Apv 2 

9 r 1 



l-2S + ¥ )K(k)-[l + 1:s )E(k) 



fc 2 



Z2pv 2 \ 3k' 2 



1 

+ W 4 



-U^ 4 - Sv 2 ) + bu 4 + 4v 2 + 8 - fc 2 (4 + 3i/ 4 ) - ^ 
Ar k 

1 



K(k) 
\2v 4 ^ 



-8fc 4 + (2 + 8i^ 2 + 4i/)fc 2 - 8 + 2v 4 + v 2 - - r r(2v 4 + I2v 2 ) + 

k l fc 4 



fc2= (1-6)(1 + C ) 



c+1 
c- 1 : 



(l + 6)(l-c)' [(l-c)(l + b)]i 



(3.1) 



(3.2a) 



E(k) , (3.2b) 



(3.3) 



Similarly the first- and third-order contributions to (3 are 



= *L[ v 2 E {k) + (fc 2 - u 2 )K(k) + (v 4 - fc 2 )n(^ 2 , fc)], 

gv 2 



and 



(3^ 



gCK(k) , «? 



l-^)x(fc)+[fc 2 (2fc 2 -l-2, 2 ) + , 4 ]|^} 



2p Apv 2 (l - c) 

_^_{_L [(2 + „ 1 + 2 „ V _ 3(;1 _ 2 „ 1KW 

^ ^ '.4 1 ,.2 ,\ 1 ;.,2/,,4 ,,2 - ~n - vA 



3fc 



+^(-(^ + ^-l) + ^(^-^+2) + -]S(fc) 



(3.4a) 



(3.4b) 



with 



v 2 = 



1-6 



1-c' 



and fc 2 and g defined in fl3.3| ). 

The first- and third-order contributions to K{= ttK) are 



;tK = -^[(fc 2 - u 2 )K(k) - v 2 E{k) - {v 4 - 2v 2 + k 2 )Il(v 2 , fc)], 



gk 2 



(3.5) 



(3.6a) 



and 



ttK2 = - 



CgK(k) 



2p 2(c 2 - l)p 

{(-2 + fc 2 + v 2 -v 4 )K{k) 



9/y 2 1 

— (fc 2 - v 2 )K(k) + j^{k 2 (k 2 -v 4 - 2v 2 ) + 2 V 4 }E{k) 



24v 4 p(c - 1) 
+2[fc 4 - (1 + v 2 + v 4 )k 2 -v 2 + 2v 4 + l]E(k)} , 



(3.6b) 



with 



and g defined in (EO). 



b- 



k 2 = 



2(6 -c) 



b+V (l + 6)(l-c) : 



(3.7) 



3.1.2. The quantities L and jj pertaining to the ^-equation [Fig. 1 in Ref. 11] 

The formulas for the first- and third-order contributions to L and L' are the same as those derived and presented 
in subsection 4.1.2 of @. 

3.2. Two complex conjugate zeros of Q 2 (rj) 
The situation of two complex conjugate transition zeros can occur only for Q 2 {rj) but not for Q 2 (£). 



3.2.1. The quantities a,f3 and K pertaining to the r\ equation: Superbarrier case 

[Fig. 4(b) in Ref. 11] 



Specializing the general formulas (2.19) and ( 2.22| ) by putting a = —1, d = +1, we obtain 



P 



1 



ijW = *iE(k)-K(k) , , 
9 l-v 2 



n 



.p 



K(k') - E(k') 



rll 



V 2 - 1' 
1 



V 2 {\-V 2 ) ^ 2 

fc 2 + (1 - k 2 )u 2 



(2j + 1)tt 



1 - v 2 \l-v 



k' 



and 



= 4- p {{-^ + ^ 1 )K(k) + ^ l (xKa ) + YEa ) 



16vj 



where 



with 



and 



" 4C + J^j) K ( k ') + 7^7 W k ') + nK(k') - E(k')}) 



A-B 1 - A-(A-B) 2 

v = -. ) H = k = —. — 

A + B y JAB' AAB 



A =[(l-h) a + <%]*, B = [(l + b 1 ) 2 +ajp 



X = -^{k 2 (U + 9v 2 ) + l7v 2 }, 



Y = —2-72 [-4A; 4 (1 + v 2 ) + fc 2 (5 + 21^ 2 ) - 17v 2 } 

On, n, 



(3.8a) 



(3.8b) 
(3.9) 

(3.10) 

(3.11a) 
(3.11b) 



We now have 



= ReffW with j = 1 



(3.12a) 
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a (3) = ReH (3) with j = 1, (3.12b) 

/3 (1) = Reff (1) with j = -1, (3.13a) 
/3 (3) = Rei? (3) with j = -1, (3.13b) 

vri? = -2IimffW, (3.14a) 

7rif 2 = -2Imff (3) . (3.14b) 
The integral 1/ associated with the contour A/y in Fig. 4(b) in |IlJ is obtained from the formula L' = a + {3. 

4. NUMERICAL ILLUSTRATION OF THE ACCURACY OF THE QUANTIZATION CONDITIONS 

For the numerical illustration of the asymmetric case we have chosen Z\ = l and considered three different values of 
Z 2 , viz. Z 2 = 2, 5 and 8. The corresponding physical systems are the ions peHe 2+ , peB 5+ and pe0 8+ , respectively, 
where p is a proton and e is an electron. For each one of these systems we have calculated the eigenvalue p and the 
reduced separation constant A' for two different c-states and various values of 7*12 ■ 

For the ion peHe 2+ we have calculated the eigenvalue p and the reduced separation constant A' for the lsa and 
2pa states and various values of ri2, with appropriate quantization conditions and parameters. The quantization 
conditions in jy] for the lsa state are (3.5a) with s = [Fig. 1(a) in Ref. 11] and (3.9) with s = m = [Fig. 4(b) in 
Ref. 11] when 7*12 is sufficiently small, but (3.5b) with s — m = [Fig. 1(b) in Ref. 11] and (3.25b) with sp = m = 
[Fig. 3(b) in Ref. 11] when 7*12 is sufficiently large. The quantization conditions in for the 2pa state are (3.5a) 
with s = [Fig. 1(a) in Ref. 11] and (3.9) with s — 1 and m = [Fig. 4(b) in Ref. 11] when t\q, is sufficiently small, 
but (3.5b) with s = m = [Fig. 1(b) in Ref. 11] and (3.25a) with s a = m = [Fig. 3(b) in Ref. 11] when r 12 is 
sufficiently large. 

For the ion peB 5+ we have computed the eigenvalue p and the reduced separation constant A' for the lsa and 
the 3scr states and various values of 7*12, with appropriate quantization conditions and parameters. The quantization 
conditions in [fllj for the lsa state are (3.5a) with s = [Fig. 1(a) in Ref. 11] and (3.9) with s = m = [Fig. 4(b) in 
Ref. 11] when r\i is sufficiently small, but (3.5b) with s = m = [Fig. 1(b) in Ref. 11] and (3.25b) with sp = m = 
[Fig. 3(b) in Ref. 11] when 7*12 is sufficiently large. The quantization conditions in [llj for the 3sa state are (3.5a) 
with s = 2 [Fig. 1(a) in Ref. 11] and (3.9) with s = m = [Fig. 4(b) in Ref. llfwhen 7-12 is sufficiently small, 
but (3.5b) with s = m = [Fig. 1(b) in Ref. 11] and (3.25b) with sp = m = [Fig. 3(b) in Ref. 11] when r X2 is 
sufficiently large. 

For the ion pe0 8+ we have calculated the eigenvalue p and the reduced separation constant A' for the lsa and 
the Ada states and various values of r±2, with appropriate quantization conditions and parameters. The quantization 
conditions in [jlTJ for the lsa state are (3.5a) with s = [Fig. 1(a) in Ref. 11] and (3.9) with s = m = [Fig. 4(b) in 
Ref. 11] when 7*12 is sufficiently small, but (3.5b) with s = m = [Fig. 1(b) in Ref. 11] and (3.25b) with sp = m = 
[Fig. 3(b) in Ref. 11] when r\% is sufficiently large. The quantization conditions in |ll| for the Ada state are (3.5a) 
with s = 1 [Fig. 1(a) in Ref. 11] and (3.9) with s = 2 and m = [Fig. 4(b) in Ref. 11] when r\i is sufficiently small, 
but (3.5b) with ~s = m = [Fig. 1(b) in Ref. 11] and (3.25b) with sp = 2 and m = [Fig. 3(b) in Ref. 11] when r 12 
is sufficiently large. 

In the calculations for the above mentioned three ions we used the quantization conditions expressed in terms 
of complete elliptic integrals, obtained in the present paper, that correspond to the above mentioned quantization 
conditions in O]. In subsection 4.1 we use the same procedure as in subsection 5.1 of |l(J for optimizing the accuracy 
of the results obtainable in the first- and third-order approximations. Thus we determine C and C as functions of ri 2 
such that the first- and third-order quantization conditions give the same results, and with these values of C and C 
we calculate p and A'. In subsection 4.2 we determine C and C such that the phase-integral quantization conditions 
give the numerically exact values of p and A' obtained by Winter et al. 0] for the ion peHe 2+ and by Ponomarev 
and Puzynina[|l3| for the ions peB 5+ and peO s+ . It is seen that the values of C and C obtained in subsection 4.1 are 
in qualitative agreement with those obtained in subsection 4.2. 
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4.1. Determination of G(r\2) and G(r\2) such that the first- and third-order quantization conditions give the 

same values of p and A' 

We have determined the values of C and C for several values of rn such that the first- and third-order quantization 
conditions give the same value of p as well as of A' . These values are tabulated and compared with the numerically 
exact results obtained by Winter et al. for the ion peHe 2+ and by Ponomarev and Puzynina|ll| for the ions peB 5+ 
and peO s+ . In Table I and Table II we give the results for the system peHe 2+ . Table III and Table IV present the 
results for the ion peB 5+ . For the ion peO s+ the results are tabulated in Table V and Table VI. The results in Tables 
I - VI are presented graphically in Figs. 1-6. 

4.2. Determination of C(r\2) and C(r\2) such that the phase-integral quantization conditions reproduce 

numerically exact values of p and A' 

By determining C and C for each value of r\2 such that the first-order quantization conditions reproduce the 
numerically exact values of p and A 1 obtained by Winter et al. Q for the lsa and 2pa states of the ion peHe 2+ and 
by Ponomarev and Puzynina|l3| for the previously mentioned two states of the ion peB 5+ and two states of the ion 
pe0 8+ we have obtained the values of C and C presented in Tables VII - XII. The numerical results in these tables 
are presented graphically in Figs. 7 - 12. 
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TABLE I: For the state Isa of the ion peHe 2+ (Zi = 1, Z2 = 2) the values of C and C have been obtained from the requirement 
that the first- and third-order phase- integral results coincide for p as well as for A' . With the use of these values of C and C 
the values of p and A' have then been obtained from the quantization conditions that are appropriate depending on whether 
ri2 is sufficiently small or sufficiently large. The numerically exact values (accurate to all digits quoted) calculated by Winter, 
Duncan and Lane [4], and obtained as private communication from Professor Winter (see p. 288-289 in [4]), are given in the 
columns called Pwdl and A' WDL . 



ri2 


C 


C 


P 


Pwdl 


P - Pwdl 


A' 


A' 

^WDL 


A' A' 


Sufficiently small ri2 


0.2 


0.4296362900 


0.4986159300 


0.2913835957 0.2909534228 


0.000430173 


-0.0401742480 


-0.0495531186 


0.00937887 


0.4 


0.4938610210 


0.4991342047 


0.5546735631 


0.5544040477 


0.000269516 


-0.16745905 


-0.1752443935 


0.007785348 


0.6 


0.5014529650 


0.5021936418 


0.7949730256 


0.7945061056 


0.00046692 


-0.3404297821 


-0.3475381522 


0.00710837 


0.8 


0.5023672195 


0.5045291831 


1.018520819 


1.018366017 


-0.000154802 


-0.54176553 


-0.5473748938 


0.00560936 


1.0 


0.5032793418 


0.5061539132 


1.231430173 


1.231534107 


-0.000103934 


-0.75725945 


-0.7624147481 


0.0051553 



2.0 0.5059145550 0.5083606740 2.241478757 2.241514227 -0.00003547 -1.870558067 -1.866548007 -0.00401006 



Sufficiently large ri2 

3.0 0.4951651150 0.5099904870 3.241389959 3.241868168 -0.000478209 -2.918225629 -2.914992386 -0.003233243 

4.0 0.4904409110 0.5094308280 4.243060447 4.243211413 -0.000150966 -3.938607817 -3.937060587 -0.00154723 

5.0 0.4878197130 0.5085272790 5.244268024 5.244326655 -0.000058631 -4.950729976 -4.949835242 -0.000894734 

6.0 0.4861586190 0.5076632010 6.245131490 6.245159553 -0.000028063 -5.958842182 -5.958257005 -0.000585177 

7.0 0.4850145200 0.5069119800 7.245774185 7.245789653 -0.000015468 -6.964658345 -6.964245417 -0.000412928 

8.0 0.4841801400 0.5062737200 8.246269561 8.246278979 -0.000009418 -7.969033663 -7.968726708 -0.000306955 

9.0 0.4835456720 0.5057325250 9.246662366 9.246668544 -0.000006178 -8.972444904 -8.972207821 -0.000237083 

10.0 0.4830475910 0.5052712520 10.24698113 10.24698542 -0.00000429 -9.975179205 -9.974990618 -0.000188587 

11.0 0.4826466130 0.5048751070 11.24724480 11.24724792 -0.00000312 -10.97741991 -10.97726635 -0.00015356 

12.0 0.4823171600 0.5045321160 12.24746643 12.24746878 -0.00000235 -11.97928963 -11.97916218 -0.00012745 

13.0 0.4820418880 0.5042327580 13.24765526 13.24765709 -0.00000183 -12.98087344 -12.98076598 -0.00010746 

14.0 0.4818086000 0.5039695200 14.24781804 14.24781949 -0.00000145 -13.98223227 -13.98214045 -0.00009182 



15.0 0.4818086000 0.5037364370 15.24795978 15.24796097 -0.00000119 -14.98341087 -14.983331513 -0.00007936 
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TABLE II: For the state 2pa of the ion peHe 2+ [Z\ = \,Zi = 2) the values of C and C have been obtained from the 
requirement that the first- and third-order phase- integral results coincide for p as well as for A' . With the use of these values 
of C and C the values of p and A' have then been obtained from the quantization conditions that are appropriate depending 
on whether m is sufficiently small or sufficiently large. The numerically exact values (accurate to all digits quoted) calculated 
by Winter, Duncan and Lane [4], and obtained as private communication from Professor Winter (see p. 288-289 in [4]), are 
given in the columns called pwdl and A' WDL . 

ri2 C C p pwdl p - pwdl A? A' WDL A' - A' WDL 



Sufficiently small r±2 

0.2 0.5294831952 0.5029173619 0.1517930283 0.1507994078 0.000993621 -2.009318372 -2.013114367 0.003795995 

0.4 0.5481736492 0.5102837463 0.3028910382 0.3062680406 -0.003377002 -2.049183021 -2.053824124 0.004641103 

0.6 0.5928374081 0.5279102827 0.468491904 0.4698372597 -0.001345355 -2.120325324 -2.125935288 0.0056099640 

0.8 0.6493296518 0.5380281837 0.642918371 0.6416038747 0.001314497 -2.226301818 -2.234147173 0.007845355 

1.0 0.6918276523 0.5400183926 0.819207383 0.8180287700 0.001178613 -2.382017327 -2.381560387 -0.00045694 

2.0 0.7291847142 0.5410643929 1.641039285 1.640235157 0.000804128 -3.839201840 -3.846791567 0.007589727 



Sufficiently large m 

3.0 0.7931298547 0.5426329856 2.304295844 2.303194434 0.00110141 -5.443294571 -5.444185235 0.000890664 

4.0 0.8368432961 0.5429173921 2.869430632 2.872046343 -0.002615711 -7.447320938 -7.448941809 0.001620871 

5.0 0.9439467950 0.5427393690 3.398416434 3.395848335 0.002568099 -9.524254856 -9.526950457 0.002695601 

6.0 0.7710812780 0.5422035320 3.913851362 3.901954918 0.011896444 -11.58670782 -11.61569923 0.02899141 

7.0 0.7080781570 0.5407494330 4.425841569 4.404864367 0.020977202 -13.63385383 -13.68696569 0.05311186 

8.0 0.6709709560 0.5389585440 4.933037475 4.909225930 0.023811545 -15.67683067 -15.73665989 0.05982922 

9.0 0.6449125300 0.5370377600 5.438128294 5.414760984 0.02336731 -17.71314008 -17.77106869 0.05792861 

10.0 0.6252417660 0.5351150510 5.942311467 5.920460066 0.021851401 -19.74285483 -19.79629610 0.05344127 

11.0 0.6098017850 0.5332640160 6.445985745 6.425795942 0.020189803 -21.76715693 -21.81592389 0.04876696 

12.0 0.5973476560 0.5315208600 6.949277104 6.930614764 0.01866234 -23.78726983 -23.83184680 0.04457697 

13.0 0.5870841030 0.5298992780 7.452239113 7.434919206 0.017319907 -25.80416211 -25.84512132 0.04095921 

14.0 0.5784759070 0.5284005350 7.954909402 7.938760040 0.016149362 -27.81854786 -27.85639810 0.03785024 



15.0 0.5711495800 0.5270195800 8.457321437 8.442196146 0.015125291 -29.83094919 -29.86611405 0.03516486 
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TABLE III: For the state Isa of the ion peB 5+ (Z\ — \,Zi =5) the values of C and C have been obtained from the requirement 
that the first- and third-order phase- integral results coincide for p as well as for A' . With the use of these values of C and C 
the values of p and A' have then been obtained from the quantization conditions that are appropriate depending on whether 
r i2 is sufficiently small or sufficiently large. The numerically exact values obtained by Ponomarev and Puzynina [13] are given 
in the columns called ppp and A' PP . 

ri2 C C p Ppp p Ppp A ^pp A App 



Sufficiently small r\2 

0.2 0.51832 0.48329 0.576736 0.57180 0.004936 -0.102639 -0.106339 0.003694 

0.4 0.51568 0.48937 1.093781 1.09319 0.000591 -0.318422 -0.312917 0.005505 

0.6 0.51379 0.49583 1.597092 1.59754 -0.000448 -0.539704 -0.533666 -0.006038 

0.8 0.51057 0.49961 2.097826 2.09827 -0.000444 -0.753875 -0.749331 -0.004544 

1.0 0.50847 0.50137 2.598201 2.59847 -0.000269 -0.962648 -0.959583 -0.003065 

2.0 0.50456 0.50273 5.099079 5.09906 +0.000019 -1.980607 -1.97997 -0.000637 



Sufficiently large ri2 

3.0 0.50320 0.50231 7.599399 7.59936 +0.000039 -2.986865 -2.98666 -0.000205 

4.0 0.50230 0.50192 10.099513 10.0995 +0.000013 -3.990163 -3.99000 -0.000163 

5.0 0.50081 0.50164 12.599680 12.5996 +0.000080 -4.991936 -4.99200 +0.000064 

6.0 0.50090 0.50141 15.099774 15.0997 +0.000074 -5.993175 -5.99333 +0.000155 

7.0 0.50082 0.50124 17.599832 17.5997 +0.000132 -6.994086 -6.99428 +0.000194 

8.0 0.50081 0.50111 20.099870 20.0998 +0.000070 -7.994774 -7.99500 +0.000226 

9.0 0.50080 0.50100 22.599897 22.5998 +0.000097 -8.995344 -8.99556 +0.000026 

10.0 0.50080 0.50090 25.099916 25.0998 +0.000116 -9.995783 -9.99563 -0.000153 

11.0 0.50070 0.50083 27.599930 27.5998 +0.000130 -10.996136 -10.9960 -0.000136 

12.0 0.50070 0.50080 30.099941 30.0998 +0.000140 -11.996448 -11.9964 -0.000048 

13.0 0.50060 0.50067 32.599950 32.5998 +0.000150 -12.996740 -12.9966 -0.000140 

14.0 0.50070 0.50060 35.099957 35.0998 +0.000157 -13.996981 -13.9969 -0.000081 



15.0 0.50060 0.50060 37.599962 37.5999 +0.000062 -14.997142 -14.9971 -0.000042 



14 



TABLE IV: For the state 3sa of the ion peB 5+ (Z\ = \,Zi = 5) the values of C and C have been obtained from the requirement 
that the first- and third-order phase- integral results coincide for p as well as for A 1 . With the use of these values of C and C 
the values of p and A' have then been obtained from the quantization conditions that are appropriate depending on whether 
r i2 is sufficiently small or sufficiently large. The numerically exact values obtained by Ponomarev and Puzynina [13] are given 
in the columns called ppp and A' PP . 

ri2 C C p Ppp p Ppp A App A App 



Sufficiently small r\i 

0.2 0.46900 0.44400 0.197728 0.196723 +0.001005 0.118904 0.077870 +0.041034 

0.4 0.43100 0.42600 0.386620 0.386538 +0.000082 0.277525 0.286919 -0.009394 

0.6 0.42620 0.42280 0.572037 0.572068 -0.000031 0.570774 0.579790 -0.009016 

0.8 0.42800 0.42280 0.754729 0.754721 +0.000008 0.916385 0.921649 -0.005264 

1.0 0.43110 0.42460 0.935242 0.935213 +0.000029 1.290323 1.29312 -0.002797 

2.0 0.44550 0.43680 1.817203 1.81720 +0.000003 3.347775 3.34845 -0.000675 



Sufficiently large ri2 

3.0 0.45680 0.44762 2.680126 2.68008 +0.000046 5.536925 5.53664 +0.000285 

4.0 0.46400 0.45600 3.533310 3.53328 +0.000030 7.777719 7.77753 +0.000189 

5.0 0.46870 0.46235 4.380874 4.38088 -0.000006 10.045943 10.0461 -0.000157 

6.0 0.47320 0.46730 5.224967 5.22494 +0.000027 12.331720 12.3314 +0.000320 

7.0 0.47608 0.47119 6.066647 6.06662 +0.000027 14.627793 14.6276 +0.000193 

8.0 0.47820 0.47430 6.906648 6.90665 -0.000002 16.931213 16.9313 -0.000087 

9.0 0.48042 0.47685 7.745476 7.74546 +0.000016 19.240602 19.2405 +0.000102 

10.0 0.48202 0.47898 8.583369 8.58335 +0.000019 21.553777 21.5537 +0.000077 

11.0 0.48337 0.48077 9.420559 9.42055 +0.000009 23.870068 23.8700 +0.000068 

12.0 0.48460 0.48230 10.257201 10.2572 +0.000001 26.188828 26.1888 +0.000028 

13.0 0.48553 0.48356 11.093399 11.0934 -0.000001 28.509499 28.5095 -0.000001 

14.0 0.48661 0.48471 11.929252 11.9292 +0.000052 30.831947 30.8318 +0.000147 



15.0 0.48720 0.48570 12.764781 12.7648 -0.000019 33.155366 33.1554 -0.000034 



15 



TABLE V: For the state lsa of the ion peO s+ [Z\ = 1, Zi = 8) the values of C and C have been obtained from the requirement 
that the first- and third-order phase- integral results coincide for p as well as for A'. With the use of these values of C and C 
the values of p and A' have then been obtained from the quantization conditions that are appropriate depending on whether 
r i2 is sufficiently small or sufficiently large. The numerically exact values obtained by Ponomarev and Puzynina [13] are given 
in the columns called ppp and A' PP . 

ri2 C C p ppp p — ppp A' A'pp A' — A'pp 



Sufficiently small ri2 

0.2 0.5019300 0.4809100 0.857612 0.855323 +0.002289 -0.137649 -0.142128 +0.004479 

0.4 0.5039400 0.4943700 1.661131 1.66144 +0.000087 -0.366556 -0.361348 -0.005208 

0.6 0.5038000 0.4989000 2.461700 2.46187 -0.000170 -0.576667 -0.573808 -0.002859 

0.8 0.5033500 0.5004300 3.261920 3.26197 -0.000050 -0.782022 -0.780407 -0.001615 

1.0 0.5028000 0.5010000 4.062060 4.06206 0.00000 -0.985331 -0.984350 -0.000981 

2.0 0.5017000 0.5012000 8.062283 8.06226 +0.000023 -1.992407 -1.99219 -0.000217 



Sufficiently large ri2 

3.0 0.5013000 0.5010000 12.062372 12.0623 +0.000072 -2.994818 -2.99479 -0.000028 

4.0 0.5008300 0.5007800 16.062383 16.0624 -0.000017 -3.996148 -3.99609 -0.000058 

5.0 0.5005000 0.5007000 20.062566 20.0624 +0.000166 -4.996555 -4.99687 +0.000315 

6.0 0.5004300 0.5005300 24.062567 24.0624 +0.000167 -5.997131 -5.99724 +0.000109 

7.0 0.5005000 0.5005000 28.062564 28.0624 +0.000164 -6.997498 -6.99764 +0.000140 

8.0 0.5004000 0.5004000 32.062561 32.0624 +0.000161 -7.997818 -7.99793 +0.000112 

9.0 0.5005000 0.5004000 36.062557 36.0624 +0.000157 -8.998017 -8.99816 +0.000145 

10.0 0.5000000 0.5004000 40.062432 40.0624 +0.000032 -9.998458 -9.99834 -0.000118 

11.0 0.5004000 0.5003000 44.062551 44.0624 +0.000151 -10.998401 -10.9985 +0.000099 

12.0 0.5000000 0.5003000 48.062446 48.0624 +0.000046 -11.998716 -11.9986 -0.000116 

13.0 0.5000000 0.5003000 52.062451 52.0625 -0.000049 -12.998804 -12.9987 -0.000104 

14.0 0.5001000 0.5003000 56.062455 56.0625 -0.000045 -13.998861 -13.9988 -0.000061 



15.0 0.4999900 0.5002100 60.062459 60.0625 -0.000041 -14.998981 -14.9989 -0.000081 
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TABLE VI: For the state Ada of the ion peO s+ {Z\ = 1, Z2 = 8) the values of C and C have been obtained from the requirement 
that the first- and third-order phase- integral results coincide for p as well as for A'. With the use of these values of C and C 
the values of p and A' have then been obtained from the quantization conditions that are appropriate depending on whether 
r i2 is sufficiently small or sufficiently large. The numerically exact values obtained by Ponomarev and Puzynina [13] are given 
in the columns called ppp and A' PP . 

ri2 C C p ppp p — ppp A' A'pp A' — A'pp 



Sufficiently small m 

0.5 0.6421800 0.5329500 0.567784 0.569166 -0.001382 -6.45764 -6.466968 +0.009328 

I. 0.6301800 0.5281010 1.597590 1.16192 -0.002161 -7.813621 -7.82220 +0.008579 
2.0 0.6223000 0.5179000 2.248680 2.24564 +0.003040 -11.571664 -11.5899 +0.018236 

Sufficiently large ri2 

3.0 0.5563400 0.5069600 3.260659 3.26096 -0.000301 -15.145897 -15.1384 -0.007497 

4.0 0.5434000 0.5090000 4.265089 4.26518 -0.000091 -18.414485 -18.4123 -0.002185 

5.0 0.5357000 0.5114000 5.265926 5.26597 -0.000044 -21.567077 -21.5660 -0.001077 

6.0 0.5303000 0.5127000 6.265543 6.26552 +0.000023 -24.662029 -24.6617 -0.000329 

7.0 0.5262000 0.5131400 7.264703 7.26468 +0.000023 -27.725910 -27.7257 -0.000210 

8.0 0.5226000 0.5131000 8.263735 8.26374 -0.000005 -30.771316 -30.7710 -0.000316 

9.0 0.5204500 0.5127600 9.262851 9.26283 +0.000021 -33.804409 -33.8043 -0.000109 

10.0 0.5186100 0.5123000 10.261986 10.2620 -0.000014 -36.829957 -36.8297 -0.000257 

II. 0.5170000 0.5118000 11.261223 11.2612 +0.000023 -39.849843 -39.8497 -0.000143 
12.0 0.5152000 0.5113000 12.260537 12.2605 +0.000037 -42.865768 -42.8656 -0.000168 
13.0 0.5142000 0.5107500 13.259930 13.2599 +0.000030 -45.878744 -45.8787 -0.000044 
14.0 0.5131800 0.5102700 14.259369 14.2594 -0.000031 -48.889604 -48.8896 -0.000004 



15.0 0.5122100 0.5098200 15.258859 15.2589 -0.000051 -51.898778 -51.8987 -0.000078 
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TABLE VII: For the state Isa of the ion peHe 2+ [Z\ — 1,Z2 = 2) the values of C and C have been obtained from the 
requirement that the first-order phase-integral results, obtained from quantization conditions that are appropriate depending 
on whether ri2 is sufficiently small or sufficiently large, coincide for p as well as for A' with the numerically exact results 
(accurate to all digits quoted) calculated by Winter, Duncan and Lane [4], and obtained as private communication from 
Professor Winter (see p. 288 - 289 in [4]), and quoted in this table as pwdl and A' WDL . 

ri2 Pwdl A' WDL C @ 



Sufficiently small ri2 


0.2 0.2909534228 


-0.0495531186 


0.51829036 


0.509271872 


0.4 0.5544040477 


-0.1752443935 


0.532092811 


0.512039821 


0.6 0.7945061056 


-0.3475381522 


0.552948218 


0.516402819 


0.8 1.018366017 


-0.5473748938 


0.562018492 


0.515629306 


1.0 1.231534107 


-0.7624147481 


0.5529481616 


0.513927361 



2.0 2.24151 -1.86655 0.5406495155 0.5122852272 



Sufficiently large r\2 



3.0 


3.241868168 


-2.914992386 


0.5148124860 


0.5120384660 


4.0 


4.243211413 


-3.937060587 


0.5060386160 


0.5106252670 


5.0 


5.244326655 


-4.949835242 


0.5018822190 


0.5092751320 


6.0 


6.245159553 


-5.958257005 


0.4994075600 


0.5081813590 


7.0 


7.245789653 


-6.964245417 


0.4977688690 


0.5072953460 


8.0 


8.246278979 


-7.968726708 


0.4965834380 


0.5065545680 


9.0 


9.246668544 


-8.972207821 


0.4957006460 


0.5059509020 


10.0 


10.24698542 


-9.974990618 


0.4950401840 


0.5054199100 


11.0 


11.24724792 


-10.97726635 


0.4943322500 


0.5050906930 


12.0 


12.24746878 


-11.97916218 


0.4940316000 


0.5045504220 


13.0 


13.24765709 


-12.98076598 


0.4934859900 


0.5044230599 


14.0 


14.24781949 


-13.98214045 


0.4933158850 


0.5041397900 


15.0 


15.24796097 


-14.983331513 


0.4931492940 


0.5037628620 
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TABLE VIII: For the state 2pa of the ion peHe 2+ {Z\ — 1,Z2 = 2) the values of C and C have been obtained from the 
requirement that the first-order phase-integral results, obtained from quantization conditions that are appropriate depending on 
whether m is sufficiently small or sufficiently large, coincide for p as well as for A' with the numerically exact results (accurate 
to all digits quoted) calculated by Winter, Duncan and Lane [4] , and obtained as private communication from Professor Winter 
(see p. 288 - 289 in [4]), and quoted in this table as pwdl and A' WDL . 

v\i pwdl A' WDL C C 



Sufficiently small ri2 

0.2 0.1507994078 -2.013114367 0.5419280238 0.50829181 

0.4 0.3062680406 -2.053824124 0.552930182 0.50920298 

0.6 0.4698372597 -2.125935288 0.572019821 0.51028392 

0.8 0.6416038747 -2.234147173 0.592038261 0.52978182 

1.0 0.8180287700 -2.381560387 0.619273017 0.542910422 

2.0 1.640235157 -3.846791567 0.632937162 0.569271625 



Sufficiently large ri2 

3.0 2.303194434 -5.444185235 0.65293027 0.551927319 

4.0 2.872046343 -7.448941809 0.6829729310 0.549321938 

5.0 3.395848335 -9.526950457 0.7689678921 0.5473760702 

6.0 3.901954918 -11.61569923 0.6760630853 0.5461458940 

7.0 4.404864367 -13.68696569 0.6153175050 0.5440702840 

8.0 4.909225930 -15.73665989 0.5787297600 0.5416277770 

9.0 5.414760984 -17.77106869 0.5582392230 0.5391808520 

10.0 5.920460066 -19.79629610 0.5461971300 0.5368472710 

11.0 6.425795942 -21.81592389 0.5382876680 0.5347039250 

12.0 6.930614764 -23.83184680 0.5325850510 0.5327246020 

13.0 7.434919206 -25.84512132 0.5281541630 0.5308665610 

14.0 7.938760040 -27.85639810 0.5245807700 0.5291925630 

15.0 8.442196146 -29.86611405 0.5216764050 0.5276965330 
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TABLE IX: For the state lsa of the ion peB 5+ (Z\ = \,Zi = 5) the values of C and C have been obtained from the requirement 
that the first-order phase-integral results, obtained from quantization conditions that are appropriate depending on whether 
ri2 is sufficiently small or sufficiently large, coincide for p as well as for A' with the numerically exact results obtained by 
Ponomarev and Puzynina [13] and quoted in this table as ppp and A' PP . 

ri2 Ppp A' pp C C 



Sufficiently small r\i 

0.2 0.57180 -0.106339 0.525368 0.501193 

0.4 1.09319 -0.312917 0.523294 0.500928 

0.6 1.59754 -0.533666 0.521226 0.500636 

0.8 2.09827 -0.749331 0.516453 0.503010 

I. 2.59847 -0.959583 0.512521 0.503770 
2.0 5.09906 -1.97997 0.505470 0.503400 

Sufficiently large m 

3.0 7.59936 -2.98666 0.503530 0.502600 

4.0 10.0995 -3.99000 0.502580 0.502110 

5.0 12.5996 -4.99200 0.502040 0.501740 

6.0 15.0997 -5.99333 0.501760 0.501410 

7.0 17.5997 -6.99428 0.501420 0.501320 

8.0 20.0998 -7.99500 0.501350 0.501030 

9.0 22.5998 -8.99556 0.501140 0.501000 

10.0 25.0998 -9.99563 0.501400 0.501300 

II. 27.5998 -10.9960 0.501220 0.501230 
12.0 30.0998 -11.9964 0.501010 0.501100 
13.0 32.5998 -12.9966 0.501000 0.501110 
14.0 35.0998 -13.9969 0.500810 0.501000 



15.0 37.5999 -14.9971 0.501000 0.500800 
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TABLE X: For the state 3sa of the ion peB 5+ [Z\ = 1, Z2 = 5) the values of C and C have been obtained from the requirement 
that the first-order phase-integral results, obtained from quantization conditions that are appropriate depending on whether 
ri2 is sufficiently small or sufficiently large, coincide for p as well as for A' with the numerically exact results obtained by 
Ponomarev and Puzynina [13] and quoted in this table as ppp and A' PP . 

ri2 Ppp A' pp C C 



Sufficiently small ri2 

0.2 0.196723 0.077870 0.428610 0.440763 

0.4 0.386538 0.286919 0.440847 0.437841 

0.6 0.572068 0.579790 0.436194 0.431006 

0.8 0.754721 0.921649 0.434223 0.428242 

I. 0.935213 1.29132 0.433012 0.426215 
2.0 1.81720 3.34845 0.446970 0.437520 

Sufficiently large m 

3.0 2.68008 5.53664 0.457260 0.448050 

4.0 3.53328 7.77753 0.464410 0.456250 

5.0 4.38088 10.0461 0.469530 0.462430 

6.0 5.22494 12.3314 0.473400 0.467340 

7.0 6.06662 14.6276 0.476400 0.471350 

8.0 6.90665 16.9313 0.478790 0.474360 

9.0 7.74546 19.2405 0.480780 0.476950 

10.0 8.58335 21.5537 0.482370 0.479130 

II. 9.42055 23.8700 0.483710 0.480810 
12.0 10.2572 26.1888 0.484920 0.482270 
13.0 11.0934 28.5095 0.485910 0.483550 
14.0 11.9292 30.8318 0.486720 0.485160 
15.0 12.7648 33.1554 0.487600 0.485520 
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TABLE XI: For the state lsa of the ion peO (Z\ = 1, Z2 = 8) the values of C and C* have been obtained from the requirement 
that the first-order phase-integral results, obtained from quantization conditions that are appropriate depending on whether 
ri2 is sufficiently small or sufficiently large, coincide for p as well as for A' with the numerically exact results obtained by 
Ponomarev and Puzynina [13] and quoted in this table as ppp and A' PP . 

ri2 ppp A'pp C C 



Sufficiently small rvi 

0.2 0.855323 -0.142128 0.494838 0.483907 

0.4 1.66144 -0.361348 0.509680 0.498740 

0.6 2.46187 -0.573808 0.506980 0.501340 

0.8 3.26197 -0.780407 0.505080 0.501930 

1.0 4.06206 -0.984350 0.504010 0.501980 

2.0 8.06226 -1.99219 0.501950 0.501460 



Sufficiently large m 

3.0 12.0623 -2.99479 0.501220 0.501160 

4.0 16.0624 -3.99609 0.501010 0.500800 

5.0 20.0624 -4.99687 0.500800 0.500700 

6.0 24.0624 -5.99724 0.500800 0.500800 

7.0 28.0624 -6.99764 0.500600 0.500700 

8.0 32.0624 -7.99793 0.500500 0.500620 

9.0 36.0624 -8.99816 0.500410 0.500600 

10.0 40.0624 -9.99834 0.500400 0.500600 

11.0 44.0624 -10.9985 0.500300 0.500500 

12.0 48.0624 -11.9986 0.500300 0.500500 

13.0 52.0625 -12.9987 0.500500 0.500300 

14.0 56.0625 -13.9988 0.500400 0.500300 

15.0 60.0625 -14.9989 0.500400 0.500200 
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TABLE XII: For the state Ada of the ion peO s+ (Zi — 1, Z2 = 8) the values of C and C have been obtained from the 
requirement that the first-order phase-integral results, obtained from quantization conditions that are appropriate depending 
on whether ri2 is sufficiently small or sufficiently large, coincide for p as well as for A' with the numerically exact results 
obtained by Ponomarev and Puzynina [13] and quoted in this table as ppp and A' PP . 

ri2 Ppp A' pp C C 



Sufficiently small ri2 

0.5 0.569166 -6.466968 0.639180 0.552970 
1.0 1.16192 -7.82220 0.628280 0.547102 
2.0 2.24564 -11.5899 0.594960 0.535600 



Sufficiently large m 

3.0 3.26096 -15.1384 0.565392 0.511310 

4.0 4.26518 -18.4123 0.546141 0.510320 

5.0 5.26597 -21.5660 0.537100 0.512080 

6.0 6.26552 -24.6617 0.531130 0.513220 

7.0 7.26468 -27.7257 0.526900 0.513530 

8.0 8.26374 -30.7710 0.523520 0.513370 

9.0 9.26283 -33.8043 0.520970 0.513030 

10.0 10.2620 -36.8297 0.519020 0.512460 

11.0 11.2612 -39.8497 0.516950 0.512100 

12.0 12.2605 -42.8656 0.515510 0.511710 

13.0 13.2599 -45.8787 0.514440 0.511010 

14.0 14.2594 -48.8896 0.513870 0.510060 



15.0 15.2589 -51.8987 0.513060 0.509620 
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FIG. 1: Plots for the lsa state of the ion peHe 2+ [Z\ = \,Zi = 2) of (a) C versus ri2, (b) C versus ri2, (c) \p — pwdl\ versus 
ri2 and (d) \A' —A' WDL \ versus ri2, when C and C are determined as functions of ri2 from the requirement that the first-order 
and the third-order phase-integral results for p as well as for A' coincide. Here p and A' are the phase-integral values obtained 
in Table I, while pwdl and A' WDL are the corresponding numerically exact values (accurate to all digits quoted) calculated by 
Winter, Duncan and Lane [4] (see p. 288-289 in [4]), and quoted in Table I. There is a break in each curve between the regions 
where the quantization conditions for sufficiently small and sufficiently large values of r i2 have been used. 
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FIG. 2: Plots for the 2pa state of the ion peHe 2+ (Z\ — 1, Z2 — 2) of (a) C versus ri2, (b) C versus ri2, (c) |p — pvkdl| versus 
ri2 and (d) |A' — A' WDL \ versus ri2, when C and (7 are determined as functions of ri2 from the requirement that the first-order 
and the third-order phase-integral results for p as well as for A' coincide. Here p and A' are the phase-integral values obtained 
in Table II, while pwdl and A' WDL are the corresponding numerically exact values (accurate to all digits quoted) calculated 
by Winter, Duncan and Lane [4] (see p. 288-289 in [4]), and quoted in Table II. There is a break in each curve between the 
regions where the quantization conditions for sufficiently small and sufficiently large values of r 12 have been used. 
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FIG. 3: Plots for the lsa state of the ion peB 5+ [Z\ = 1, Z2 = 5) of (a) C versus ri2, (b) C versus ri2, (c) \p — ppp\ versus r\2 
and (d) \A' — A' PP \ versus ri2, when C and C are determined as functions of m from the requirement that the first-order and 
the third-order phase-integral results for p as well as for A 1 coincide. Here p and A' are the phase-integral values obtained in 
Table III, while p pp and A' PP are the corresponding numerically exact values obtained by Ponomarev and Puzynina [13] and 
quoted Table III. There is a break in each curve between the regions where the quantization conditions for sufficiently small 
and sufficiently large values of r i2 have been used. 
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FIG. 4: Plots for the 3sa state of the ion peB 5+ (Z\ = 1, Z2 = 5) of (a) C versus ri2, (b) C versus ri2, (c) \p — ppp\ versus 
r\2 and (d) \A' — A' PP \ versus ri2, when C and C are determined as functions of m from the requirement that the first-order 
and the third-order phase-integral results for p as well as for A' coincide. Here p and A' are the phase-integral values obtained 
in Table IV, while ppp and A' PP are the corresponding numerically exact values obtained by Ponomarev and Puzynina [13] 
and quoted in Table IV. There is a break in each curve between the regions where the quantization conditions for sufficiently 
small and sufficiently large values of rvi have been used. 
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FIG. 5: Plots for the Isa state of the ion peO s+ (Z\ = 1, Z2 = 8) of (a) C versus ri2, (b) C versus m, (c) \p — ppp\ versus 
ri2 and (d) \A' — A' PP \ versus ri2, when C and C are determined as functions of m from the requirement that the first-order 
and the third-order phase-integral results for p as well as for A' coincide. Here p and A' are the phase-integral values obtained 
in Table V, while ppp and A' PP are the corresponding numerically exact values obtained by Ponomarev and Puzynina [13] and 
quoted in Table V. There is a break in each curve between the regions where the quantization conditions for sufficiently small 
and sufficiently large values of r\i have been used. 
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FIG. 6: Plots for the Ada state of the ion peO s+ (Z\ = 1, Z2 = 8) of (a) C versus m, (b) C versus m, (c) \p — ppp\ versus 
ri2 and (d) \A' — A' PP \ versus ri2, when C and C are determined as functions of m from the requirement that the first-order 
and the third-order phase-integral results for p as well as for A' coincide. Here p and A' are the phase-integral values obtained 
in Table VI, while p pp and A' PP are the corresponding numerically exact values obtained by Ponomarev and Puzynina [13] 
and quoted in Table VI. There is a break in each curve between the regions where the quantization conditions for sufficiently 
small and sufficiently large values of m have been used. 
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FIG. 7: Plots for the lsa state of the ion peHe 2+ [Z\ = 1, Z2 = 2) of (a) C versus r\i and (b) C versus ri2, when C and C are 
determined as functions of ri2 from the requirement that the first-order phase- integral results and the numerically exact results 
obtained by Winter el al. [4] coincide. There is a break in each curve between the regions where the quantization conditions 
for sufficiently small and sufficiently large values of ri2 have been used. 
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FIG. 8: Plots for the 2pa state of the ion peHe 2+ (Zi = 1, Z2 = 2) of (a) C versus ri2 and (b) C versus ri2, when C and C arc 
determined as functions of ri2 from the requirement that the first-order phase-integral results and the numerically exact results 
obtained by Winter et al. [4] coincide. There is a break in each curve between the regions where the quantization conditions 
for sufficiently small and sufficiently large values of ri2 have been used. 
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FIG. 9: Plots for the lsa state of the ion peB 5+ [Z\ = 1, Z2 = 5) of (a) C versus r\i and (b) C versus ri2, when C and C are 
determined as functions of ri2 from the requirement that the first-order phase- integral results and the numerically exact results 
obtained by Ponomarev and Puzynina [13] coincide. There is a break in each curve between the regions where the quantization 
conditions for sufficiently small and sufficiently large values of ri2 have been used. 
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FIG. 10: Plots for the 3sa state of the ion peB 5+ (Z\ = 1, Z2 = 5) of (a) C versus ri2 and (b) C versus ri2, when C and C are 
determined as functions of ri2 from the requirement that the first-order phase-integral results and the numerically exact results 
obtained by Ponomarev and Puzynina [13] coincide. There is a break in each curve between the regions where the quantization 
conditions for sufficiently small and sufficiently large values of ri2 have been used. 
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FIG. 11: Plots for the Isa state of the ion peO a+ (Z\ = 1, Z2 = 8) of (a) C versus r\2 and (b) C versus ri2, when C and C are 
determined as functions of r\2 from the requirement that the first-order phase- integral results and the numerically exact results 
obtained by Ponomarev and Puzynina [13] coincide. There is a break in each curve between the regions where the quantization 
conditions for sufficiently small and sufficiently large values of r\2 have been used. 
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FIG. 12: Plots for the Ada state of the ion peO s+ (Zi = 1, Z2 = 8) of (a) C versus ri2 and (b) C versus ri2, when C and C arc 
determined as functions of ri2 from the requirement that the first-order phase-integral results and the numerically exact results 
obtained by Ponomarev and Puzynina [13] coincide. There is a break in each curve between the regions where the quantization 
conditions for sufficiently small and sufficiently large values of ri2 have been used. 



